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A calculation of the semi-leptonic decays of the kaon (Kis) is presented. The results are direct 
predictions of a covariant model of the pion and kaon introduced earlier by Ito, Buck, Gross. The 
weak form factors for K13 are predicted with absolutely no parameter adjustments of the model. 
We obtained for the form factor parameters: f-(q 2 = mf)/f + (q 2 = mj) = —0.28 and A + = 0.028, 
both within experimental error bars. Connections of this approach to heavy quark symmetry will 
also be discussed. 
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I. INTRODUCTION 



The success of the Ito-Buck-Gross (IBG) ]lj model in the description of many properties of both the ir and K 
mesons motivated the calculation of the Ki 3 decays reported on here. It is the Ki 3 decays that combine both the pion 
■ and kaon wave functions generated previously [1-3]. A successful K13 calculation that is coupled to other observables, 
\Q ' constrains further the physics described by the model. 

The work reported on here is predictive and employs no new parameters and no parameter adjustment. The results, 
found below, are in good agreement with the data and are as good as Chiral Perturbation Theory (CPT) approaches 
H and ||, even though a low energy parameter is involved in them, and better, at least in the light-quark sector, 
r"| ■ than the Quark Potential Model (ISGW2), with a hyperfine interaction, predictions ||. It is noted that an older 
Qh! version of the Quark Potential Model without hyperfine interaction (ISGW) gives results similar to ours Q|. These 
comparisons and the details of our calculation are presented below. Nonetheless, with the success of CPT and this 
' work, the question still remains for nuclear physics as to singling out quark structure from hadronic. That is, where 
do hadrons leave off and quarks begin? The Quark Model has been very successful at reproducing hadronic static 
properties such as the mass spectrum and moments. But it is the dynamic properties, we feel, that will delineate 
the differences between hadronic physics and quark physics. For this reason, we take the position that not only is 
the near, or at q 2 = physics important but the q 2 ^ domain will delineate theoretical approaches. Thus, our 
predictions for the non-perturbative weak transition form factors as a function of q 2 are also presented here in an 
attempt to draw both theoretical and empirical interests. 

A detailed review of the theoretical and experimental status of semileptonic kaon decays is given in Ref . M . 



II. THE MODEL 



The theoretical model employed is an extension of the Nambu-Jona-Lasinio (NJL) model but with a definite 
momentum distribution generated by a separable interaction: 

V(p, k) = gf(p 2 )f{k 2 )[I (81- ( 7 5 r) ® ( 7 5 r)], (1) 
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where f(p 2 ) — (A 2 — p 2 )^ 1 with A being the interaction cutoff parameter for a given meson state. With this choice 
of the qq interaction, one can integrate all momentum integrals to infinity; there is no need for an integral cutoff as 
employed in the NJL model. 

The IBG model requires that the Bethe-Salpeter equation be solved for the vertex function, T, for each meson 
considered, 



P 




k-p/2 



where p is the 4-momentum of the meson. 

The self energy of each flavour quark is treated by solving the Schwinger-Dyson equation 

where n/ is the number of quark flavours (equal to 3 in our qq system) in a coupled sense (coupled via quark masses 
and interaction strength) to the Bethe-Salpeter equation. Though, in the case of the strange quark mass, the self 
energy is assumed to be the (constituent) quark mass and is treated as a parameter to be fixed. |J 
Electomagnetic gauge invariance is imposed on the electromagnetic current of a pseudoscalar meson 

J» = F(q 2 )(p + p'r, (3) 

where F(q 2 ) is the meson charge form factor and p(p') is the 4-momentum of initial (final) meson. The work of Buck, 
Williams, and fto (BWI) || has shown that both the pion and kaon (charged and neutral) charge form factors can be 
predicted and it is the pion and kaon vertices from this work that are employed in the calculation of the weak form 
factors. 



III. WEAK FORM FACTORS 



In the Standard Model, the weak current for K13 decays has the following structure 

= ^v us {f + ( q 2 )(p K + p„r + f-(q 2 )(PK - p*Y] (4) 

where Pk and P T are the kaon and pion 4-momenta, q — Pk — P-n, and f± are dimensionless form factors. 
The semi-leptonic decays studied are: 

{K&)K±=-n a e k v e 
K° L = n ± eTv e 

(K, a )K ± =Ti ii ± u^ (5) 

Kl = T^/X^ 

In the limit of exact isospin symmetry, m u = m c i, form factors of charged and neutral kaon decays are related: 

fi/fi = 1/V2, 

and in the limit of exact SU(3) symmetry, the form factor /_ is zero. For the decay channel, the transferred 4- 
momentum q is time-like, and the physical region is limited to mf < q 2 < (txik — m-n) 2 . The vertices appearing in this 
weak current and the ones employed in this work are the kaon and pion vertices (wave functions) previously obtained 
by BWI; namely 
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r ^-( fc ) = A | g '!_ 7 fc2 » ( 6 ) 



with A^if j7r being the normalization. 
From equation 4, one can uncouple f±: 



{Pk T P^) 2 J"(Pk ± P*)" - (m^_- ml)J»(P K T P*Y 

■( 

To compare to available experimental data, the following low-g 2 expansion is used for the form factors: 



J± ~ (Pk - P,) 2 (Pk + P.) 2 - Wk - ™IY 



/±(9 2 ) = f±(q 2 = mf){l + A ± ^^) (8) 

m t+ 

where A± is the slope of f± evaluated at q 2 = mf and f±(q 2 = mf) corresponds to the normalization. Note that it is 
the charged, not the neutral, pion mass that enters the above expansion. 

Another set of the form factor parameters commonly used in the literature is A + , Ao, arising as coefficients of linear 
expansion of the form factors /+ and /rj, with fo defined as 

/o = /++ 2 q 2 /-, (9) 

The form factors /+ and fo describe, respectively, P— wave and S— wave projections of weak current matrix elements 
in the crossed channel. 

To obtain the values of A± , a calculation of J M must be performed. In this work, J M is the direct result of a triangle 
diagram (Fig.l) with a flavour changing operator having V — A spin structure 7^(1 — 7 5 ). In the Standard Model, 
the K13 decay form factors are determined only by the vector part of the charged weak current operator. 

W(q) 



^ = 



K(P K ) 



FIG. 1. Triangle diagram for the charged weak current of K13 decay. 

Integrals with respect to loop momentum were evaluated in the following way: In the expession for the weak current 
given by a Feynman diagram Fig.l, the spin trace was calculated and the terms dependent on loop momentum in 
the numerator were divided out by corresponding terms in the denominator. This procedure reduces the expression 
for 'impulse' current Fig.l to the sum of scalar integrals of products of three to five denominator factors (three of 
them coming from quark propagators and two from meson-gg vertex form factors). Each denominator factor is a 
polynomial quadratic in the loop momentum. The terms involving three denominators, are in fact scalar 3-point 
functions which may be expressed analytically in terms of Spence functions In this work, to calculate the 3-point 
functions, we parametrized them in terms of two Feynman parameters. Integration with respect to one Feynman 
parameter was done analytically, and the other was a numerical Gauss integration. We did not use a low q 2 expansion 
to evaluate loop integrals but we do extract from our results the low q 2 behavior and find it is consistent with the low 
q 2 expansion of equation 8 above, employed by all researchers. 

Our task of computing 4-dimensional loop integrals of products of more than three denominators is greatly simplified 
when taking advantage of the fact that only two external momenta in the integrand are linearly independent. As a 
result, products of four and five denominators are reduced to the sum of products of three denominators with redefined 
masses Mi. This procedure is described by the following identity 

' (10) 



(k 2 - m 2 )((k + qi) 2 - m 2 )((fc + q 2 ) 2 - m 2 ){{k + q 3 ) 2 - m 2 ){{k + q 4 ) 2 - m 2 ) 

E 



i,j.n 



(k 2 - M 2 )((k + q 3 ) 2 ~ M 2 )((k + q n ) 2 - Ml) 
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if 

(:0= A (£)' detA *°' 

where the sum is taken over different combinations of the external 4-momenta qi involved in the reaction, m, are 
quark masses and mass parameters in meson-gg vertex form factors, a.; are coefficients independent of loop momentum 
k, and A is a 2x2 matrix setting relations between external momenta in the integrand. After this reduction, scalar 
3-point integrals are computed within the technique described above. 

IV. RESULTS 

In the physical region of K e3 decays, q 2 may be as low as the square of electron mass, 25-10 -8 GeV 2 /c 2 , and as 
high as the square of the mass difference between the kaon and the pion, 0.123 GeV 2 /c 2 . The form factors f± in this 
region with a good precision appear to be linear functions of q 2 thereby justifying a linear parametrization of equation 
8 usually employed in analyses of experimental data JlOj] . To compare our results with experiment, we extracted the 
slopes and ratios of the form factors f± at q 2 = mf via numerical differentiation. Numerical values for the parameters 
in this calculation were taken to be the same as in Ref. ||, viz. m Ui d— 250 MeV, m s = 430 MeV, A^— 600 MeV, and 
A K = 690 MeV. 

The direct predictions of our approach for A+ and j^\ q 2 =m 2 are 0.028 and —0.28, respectively. These results are to 
be compared to the experimental values of 

A + = 0.0286 ± 0.0022 and U =/_//+ = -0.35 ±0.15. 

We obtain A_= 0.029, i.e., in our model both /_ and / + have approximately the same slopes, in agreement with early 
Quark Model results |T^ |. Our calculation for K e a and yields equal results, within the quoted precision, since 
the A± are almost constant in the range m 2 < q 2 < m 2 . Naturally, the decay rates should be different due to phase 
space factors; they can be calculated by known formulas in terms for form factor slopes, see, e.g. Ref. 0; however, 
we have yet to perform the calculation of these rates. 

Table 1 illustrates the comparison between our work, that of CPT, vector meson dominance (VMD), and the ISGW2 
model. One sees that the work reported on here compares very favorably to experiment and to CPT, except for the 
prediction for Ao- 

Tablel: Model predictions for the parameters of decay form factors. 

* From the corresponding values of A+ and Aq 0; ** From the corresponding values of A + and £a(0). 





CPT: Refs. §] and 
p|, respectively 


VMD §] 


ISGW2 @ 


This work 


Experiment 


A+ 


0.031; 0.0328 


0.0245 


0.019 


0.028 


0.0286±0.0022 (K e3 ) 


6t(o) 


-0.164±0.047*;-0.235 


-0.28 


-0.28 


-0.28 


-0.35±0.15 (Jfpa) 


Ao 


0.017±0.004; 0.0128 


0.0 


-0.005** 


0.0026 


0.004±0.007 (K+ 3 ) 
0.025±0.006 (JT° 3 ) 



A prediction for the slope parameter Ao obtained within our model is 0.0026, which is consistent with experiments on 
charged kaon decays (Aq = 0.004±0.007) and inconsistent with neutral kaon decay measurements ( Aq = 0.025±0.006). 
Since the experimental results for this slope parameter are not firm, it is hard to draw any positive conclusions about 
agreement or disagreement of our result for Ao with experiment. However, we can compare it with predictions from 
other models. It can be seen from the Table 1 that the quark model in general gives much smaller numbers for Ao 
than CPT. 

To test the sensitivity of our results, an arbitrary change in the A^ cutoff from 600 MeV to 450 MeV (a 25% 
change), results in a A + = 0.028 to 0.031 (an almost 10% change), £a(0) from —0.28 to —0.27 (« 5% change), and 
Ao from 0.0026 to 0.0058 (« 100% change), respectively. One is reminded that changing A^ changes the pion charge 
radius as well as the pion decay constant. In fact, the value A w = 450 MeV was used in Refs. [1,2] as the best fit to 
pion decay constant and charge radius alone, and this parameter was adjusted to 600 MeV in Ref. H to be able to 
treat both pion and kaon in a coupled approach. 

It should be noted that our model gives a stable, with respect to variation of the model parameter, prediction for 
£a(0), and appears to give a highly parameter-dependent result for Ao. This model dependence is due to cancellation 
between two large terms on the r.h.s. of equation (9). The situation is different in CPT, where uncertainties due to 
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higher-order loop corrections give rise to about 30% uncertainty for £a(0): — 0.164±0.047 and about 25% uncertainty 
forAoi- 

Finally, CPT suggests that once the ratio of the weak form factors is known, then an estimate of the mass of the 
strange sigma (m aK ), a meson with J p = + , can be made. The relationship referred to is H 



U(0) = {M 2 K - Ml){M p 



M: 



(ii) 



Taking our result for £a(0) and assuming duality between our model predictions and the model with effective 
exchanges of vector and scalar mesons at low q 2 , we have m aK — 1.5 GeV, which compares favorably to the mass 
of -?Cq(1430). A test of this value could be made through the hypernuclear spectroscopy measurements (CEBAF 
E89-009; CEBAF PR-95-002) inferring the interaction that contains this strange sigma. 

Another feature of our approach is revealed in the limit as A^, m^-, and m s become infinitely large. The ratio 
d(0) is calculated in this limit and its asymptotic value is -1; Figure 2 illustrates the mass dependence. In the limit 
of Heavy Quark Symmetry (HQS), the q 2 dependence of the semileptonic decay form factors is factorized out in the 
form of Isgur-Wise function |11|, and £a is given by the combination of the initial (M) and final (M') meson masses, 
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M-M' 
M + M' ' 



(12) 



As a result, £,a\hqs = ~ 1 if the initial meson is much heavier than the final. Note that for HQS to be applied, both 
initial and final mesons should be heavy, whereas assuming m s to be large in our model, we keep the final meson 
light. This implies that this particular result of HQS appears to be more generally applicable. Of course, the ratio is 
zero for mesons of equal mass. 




10.0 
(M-m K ),GeV 



FIG. 2. Ratio £a(0) as a function of initial meson mass. 



Though, it is tempting to make exuberant statements in regard to the identical results, one is cautioned by the 
manner in which the limits are taken, and to the nature of the physics examined, respectively. 

Furthermore, Figure 3 illuminates our predictions for f + (q 2 ) at space-like momentum transfers describing the 
neutrinoproduction processes vir —* IK, vir — > lir + and corresponding weak lepton capture. We stress that no low 
q 2 expansion was assumed in our calculations, so that present results have the same validity range in terms of q 2 as 
results of Refs.[l-3] for electromagnetic form factors of pion and kaon. The form factor / + at large q 2 behaves as 
1/q 4 (up to a logarithmic correction) indicating that our model effectively describes soft, nonperturbative reaction 
mechanism, and does not include perturbative QCD contributions. 

The weak K — > 7r transition form factors in the space-like region could be possibly accessed experimentally in the 
production of kaons on an hadronic target induced by neutrinos or lepton weak capture. The latter possibility is being 
studied for a CEBAF experiment Q. 

It would be instructive to see if the earlier success of the IBG model, that include the pion and kaon observables as 
well as the results of this present work, can be reproduced with other interactions and/or with other wave equations; 
by this, it is meant the predictive characteristics associated with the low energy axial anomaly, such as the pion 
transition and elastic charge form factors, the kaon charge form factors, and the decays. 

In brief summary, weak form factors and slope parameters have been calculated for decays. The results compare 
very favourably to available experimental data. The model employed was that of IBG Refs.[l-3] and there were no 
parameter adjustments, thus, rendering this calculation predictive. 
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